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We argue that a generic instability afflicts vacua that arise in theories whose mod- 
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fields we provide numerical evidence that for a generic local minimum of the poten- 
tial the usual semiclassical bubble nucleation rate, F = Ae~^, increases rapidly as 
function of the number of fields in the theory. As a consequence, the fraction of vacua 
with tunneling rates low enough to maintain metastability appears to fall exponen- 
tially as a function of the moduli space dimension. We discuss possible implications 
for the landscape of string theory. Notably, if our results prove applicable to string 
theory, the landscape of metastable vacua may not contain sufficient diversity to 
offer a natural explanation of dark energy. 
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I. INTRODUCTION 

The discovery that string theory admits an enormous number of flux vacua [liJal has 
played a formative role in the theory's development for well over a decade. Aspects of 
these vacua, from their phenomenological and cosmological properties to their distribution 
and statistical features, have been extensively studied Q, Isj. Yet, due to the complexity of 
this landscape of vacua, many basic questions remain. In this paper we undertake a study, 
pursued from a variety of perspectives in a number of works (6|-[l3j, that is of potential 
relevance to one such vital question: Do we expect these vacua to be long-lived? As a direct 
analysis would present formidable challenges, we instead consider generic, field theoretic 
models of the landscape and study how the stability of vacua varies as the dimension of the 
moduli space (the number of fields) increases. Our results suggest that tunneling rates, and 
hence vacuum instability, grow so rapidly with the number of moduli that the probability 
of a given local minimum being metastable is exponentially small. 

In field theory, vacuum decay by quantum tunneling was studied by Coleman HQ. 
He showed that the decay proceeded by the nucleation of bubbles of a lower vacuum inside 
the original false vacuum. In the semiclassical approximation the nucleation rate per unit 
volume is governed by a bounce solution of the Euclidean field equations. It can be written in 
the form F = Ae~^ , where A depends on the determinant of fluctuations around the bounce 
solution and B is the Euclidean action of the bounce. The analysis was extended to include 
gravitational effects by Coleman and De Luccia [16j. For decay from a de Sitter vacuum 
the resulting corrections to the nucleation rate are typically small unless the potentials are 
Planckian in scale or the bubbles nucleate with a size comparable to the horizon length. 
With unusually flat potential barriers it can happen that there is no Coleman-De Luccia 
bounce, but in such cases there is always a Hawking-Moss 17|] solution corresponding to a 
process in which an entire horizon volume fluctuates to the top of the potential barrier. 

In string theory the large number of vacua and moduli flelds complicates the situation, 
but at the same time opens up lines of attack based on statistical analysis. For example, 
Denef and Douglas 18| proposed a method of calculating the density of flux vacua in the 
string landscape in terms of the Kahler potential on the moduli space of a given Calabi-Yau 
compactiflcation. Their work showed that a sharp accumulation of vacua generally occurs 
near the conifold locus in moduli space. Dine et al. jsl used scaling arguments to conclude 
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that vacua in the string landscape with small cosmological constant become unstable when 
fluxes are large relative to their compactiflcation volume. Chen et al. [lol and Marsh et 
al. 20[ showed that increasing the number of moduli fields suppresses exponentially the 
chance that a randomly chosen critical point will be a minima, with the suppression growing 
for vacua with high energy. 

In this paper we also pursue a statistical approach, focusing our analysis on effective field 
theory models of many- dimensional moduli spaces. We provide numerical evidence that the 
rate for tunneling out of a typical false vacuum grows rapidly as a function of the number of 
moduli fields. Specifically, the fraction of vacua with tunneling rates low enough to maintain 
metastability appears to fall as an exponential of a power of the moduli space dimension. 

In Sec. miwe describe our approach for estimating tunneling rates in field theoretical mod- 
els of high- dimensional moduli spaces. Our numerical methods and results are described in 
Sec. mil These results reveal a general feature of high-dimensional field theories and are in- 
dependent of applications to the string landscape. In Sec. [IV]we discuss the efficacy of using 
random potentials to model the string landscape, and emphasize various considerations that 
would need to be resolved before a direct application could be justified. In Sec. |V]we esti- 
mate the maximum dimension of moduli spaces whose associated fiux landscape would be 
expected to generically contain metastable vacua, assuming that our numerics and extrapo- 
lations are applicable. The result is a drastic reduction in the number of metastable vacua. 
Finally, in Sec. I VI I we suggest future directions for studying these models and summarize 
our conclusions. 



II. VACUUM DECAY 



We consider the dynamics of a moduli space spanned by N scalar fields with a La- 
grangian 

1 ^ 

^ = 2 E - V{<Pu 02, ... , <Pn) . (2.1) 

i=i 

The potential V will in general have many local minima that correspond to metastable false 
vacua. Let us consider one of these which, by a shifting of the field variables, can be taken 
to lie at the origin of field space, = 0. Assuming the potential to be smooth at this point, 
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we can expand it in a power series 

^ = ^ ( E + E ^l^i^^^'^^ + E + - ) ■ (2.2) 

\ i ijk ijkl / 

Here with dimensions of mass, is a characteristic distance in field space that corresponds 
to a typical distance between stationary points of V . We have also extracted a dimensionless 
constant A, to be chosen so that the dimensionless coefficients of the power series in brackets 
are of order unity. Finally, we have used the freedom to make an 0{N) transformation on 
the fields to eliminate off-diagonal terms in the quadratic term of the power series. 

The exponent B in the bubble nucleation rate is the action of the Euchdean bounce 
solution. We ignore gravitational effects and assume 0(4) symmetry, with the fields being 
functions only of s = a/x^ + xl- The bounce then satisfies 

ds^ s ds d(j)j 

The boundary conditions are that 0(oo) = 0, its false vacuum value, and that 0'(O) = 0. 
The actual value of the field at the origin is not determined in advance, but must be a point 
on the opposite side of the potential barrier from the false vacuum. Note that, except in the 
thin-wall limit, 0(0) is never equal to the true vacuum value. 

We will consider large ensembles of potentials, with the coefficients in the power series 
chosen randomly, as described in more detail in the next section. Ideally, we would calculate 
the nucleation rate for each potential by solving the bounce equations. However, finding 
bounce solutions in a theory with more than one scalar field is a daunting numerical problem. 
Doing so for a large sample of potentials is clearly infeasible. Instead, we invoke a more easily 
calculable proxy that can provide an indication of how the decay rate varies with the number 
of fields. 

A bounce solution may be viewed as containing a wall region, in which the fields pass 
through the barrier in T^(0), and an interior region where the fields are close to their values 
in the new vacuum. The distinction between the two regions becomes exact in the thin-wall 
approximation, which is valid in the limit where e, the difference between the energy densities 
of the true and false vacua, tends to zero. In this approximation the fields in the interior 
region take on exactly their true vacuum values. The radius R of the bounce is determined 
by a balance between the negative action in the bounce interior and the positive action in 
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the wall. With a single scalar field, this wall action is the product of the three-dimensional 
area of bounce and a surface tension that is well approximated by 



where (j)^, denotes the point on the true vacuum side of the barrier such that V^cj)^) = V{^. 
With additional scalar fields, a is obtained by integrating along a path in field space running 
from the false vacuum to a point near the true vacuum, with the path and endpoint chosen 
to minimize the integral. 

The net result is that R = 3a/e, while the tunneling exponent is 



Thus, in the rare cases in which we have two almost degenerate vacua, the thin-wall ap- 
proximation is applicable, B is large, bubble nucleation is greatly suppressed, and the false 
vacuum is long-lived. Since our focus is on effects that enhance bubble nucleation, this 
approximation is not of direct interest to us. Nevertheless, we can draw some useful insight 
from it. Equation (12. 5 p shows a strong dependence on e. This cannot continue outside the 
thin-wall approximation, because then the field in the bounce solution never reaches the 
true vacuum, and so cannot be directly sensitive to the value of e. On the other hand, the 
surface tension is closely related to the form of the potential barrier and should continue to 
be relevant. 

Outside the thin- wall limit, the boundary between the wall and the bounce interior is 
not well defined. A reasonable prescription would be to take it to be the hypersurface in 
Euclidean space on which ¥{(/)) is equal to its false vacuum value. We could then define a 
quantity a as before, with the integration path running from the false vacuum to a point on 
the hypersurface S in field space, lying on the other side of the potential barrier, on which 
V{(f)) = Vfv. The path and the specific endpoint on S would be chosen to minimize the 
integral. 

Unfortunately, performing the required minimization for a large ensemble of potentials 
is still calculationally infeasible. However, a plausible approximation is at hand. We might 
expect the minimizing path to pass though the region in field space where the barrier in V{(f)) 
is lowest. This suggests considering a straight-line path running from = 0, though each 
saddle point, 0sp, on the surrounding barrier, and ending on S. In fact, since on average the 




(2.4) 



B = -aR' 
2 



277rV 
2e3 



(2.5) 
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contributions from the segments before and after the saddle point will be equal, we follow a 
slightly simpler approach. We integrate over only the first part of the path, extract factors 
of A and f , and define the line integral 

a = d(j) ^/2[V{(P) -V{0)] 

= y/xv^s, (2.6) 

where P is a straight-line path in field space running from the false vacuum at = to the 
saddle point. There may be many saddle points, and thus many such paths. We expect the 
tunneling rate to be controlled by the path with the lowest tunneling integral, and define 
the corresponding value of s to be 

^min "5 ■ (^•'^) 

Standard scaling arguments show that B is inversely proportional to A, but independent 
of V, although v does affect the prefactor A in F. Outside the thin-wall limit, the typical 
bounce radius is i? ~ k{\/Xv)~^ , where A; is a numerical factor of order unity. This suggests 
that 

B ^ Ti'^R^d ^ —k^ s , (2.8) 
A 

where we have included a factor of tt^ because of the four- dimensional spherical symmetry. 
A simple test of this estimate can be obtained by numerically evaluating B and s for a 
single-field potential with cubic and quartic terms. Inserting the results into Eq. (12. 8p and 
taking v to be half the difference between the true and false vacuum values of 0, one obtains 
values of k that vary between 5.2 and 6.0 over a wide range of parameters away from the 
thin- wall limit. This suggests 

5-10=^^. (2.9) 
A 

The dilute-gas approximation that underlies the semiclassical approach to bubble nucle- 
ation breaks down when B is comparable to or less than unity. In this regime the metastabil- 
ity of the false vacuum has essentially disappeared. The numerical studies that we describe 
in the next section show that when many scalar fields are present, the overwhelming major- 
ity of potentials lead to a value of B too small to maintain metastability with any plausible 
value of A. 
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III. NUMERICAL STUDIES 

We numerically studied ensembles of theories with potentials of the form of Eq. (12.21) . 
using s as an indicator of the vacuum stability. To make the calculations manageable we 
truncated the power series at the quartic terms. An ensemble was defined by taking the 
A'-"^ to be random numbers uniformly distributed over ranges defined by 

4? e [0,a2], 

e [-03,03], 

4L e [-a„a,]. (3.1) 

Allowing the A^^^l^ to be negative means that the truncated potential is not necessarily 
bounded from below. This is not a concern for us, since we are only concerned with the 
behavior near the minimum; at larger distances higher-order terms can provide a lower 
bound on V. 

Because we will be comparing theories with different numbers of fields, an important issue 
is how to vary the a„ as the number of fields is varied. To determine this, let us require that 
the typical variation of the potential in a ball of radius (ppt in field space be independent of 
N, the number of fields. This ensures that, for any number A^^o < A^, we can recover results 
for N — Nq fields by considering an A^g-dimensional cross section of the analysis for A^ fields. 
In turn, this ensures that dependencies we find on A^ are not due to peculiar A^-dependent 
normalizations in V. 

With this assumption, the typical value for each of the 0j in such a ball of radius 
is of order (pj^/y/N. The quadratic term is a sum of positive contributions and will be 
independent of A^ if 02 is. There are A^^ cubic terms, each of magnitude 0|./A^^/^. Because 
these can be of either sign, they will tend to cancel, so that the effective number of terms 
is of order N^/'^ and we are led to take 03 to also be A^- independent. A similar argument 
shows that 04 should also be independent of A^. 

As for the actual values of the a„, note that a change of any two of these can be absorbed 
by a redefinition of A and v. Hence, there is no loss of generality in taking 02 = 03 = 1. We 
did so, and also set 04 = 1; the effect of other choices for 04 is described below. 

For a given value of A^ we chose an ensemble of 10,000 potentials. For each potential we 
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FIG. 1: Left, the median value of s for quartic potentials. Right, the median value of the lowest 
saddle point, in units of Af^, for quartic potentials. In both cases the bars indicate the range from 
the 25th to the 75th percentile. 

found all of the stationary points^ and picked out the saddle points with a single negative 
mode. From among the saddle points for each potential we found the one that gave the 
smallest value of s (i.e., s^in = s) and the one that gave the lowest barrier height. Figured] 
shows the median values of these quantities within each ensemble as functions of for the 
quartic potential. For both quantities a sharp decrease with increasing is clearly evident. 
Not only do the height and surface tension of these saddle points decrease, but also they 
get closer to the false vacuum minimum at the origin. This can be seen in Fig. [2], where 
we have plotted the median distance to the saddle point with lowest s. (A plot of median 
distance to the lowest saddle point is virtually indistinguishable.) All of these plots show a 
power law dependence on A^ with, e.g.. 



■^median ~ C*i 



tension 



(3.2) 



The best fit values for the various a and C are shown in Table [B 

We have also plotted the number of stationary points around the false vacuum, in Fig. [3l 
These do not follow a power law, but are instead closer to an exponential behavior. 

Our decision to arbitrarily terminate the expansion of the potential with the quartic terms 
was motivated by considerations of calculational practicality. As a test of this choice, we 
also carried out the calculations without the quartic terms in the potential, retaining only 



^ For a description of the numerical method for finding all such points, see 21 1. Further details will be 
published elsewhere Q. 
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FIG. 2: The median distance, in units of v, to the saddle point with Smin = s for quartic potentials. 
Again, the bars indicate the range from the 25th to the 75th percentile. 
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TABLE I: Best fit parameters, defined as in Eq. ()3.2p . for power law fits to the data in Figs. [H [21 
m and [5l as well as to the data for the super symmetric potentials discussed in Sec. llVi 

the quadratic and cubic terms. As can be seen from the plots in Figs. HHHl and the data in 
Table [U the results are quite similar to those with the quartic term included. This leads us 
to conclude that our omission of quintic and higher terms has little effect on our results. 

We now return to the question of the dependence of the results on the ranges chosen for 
the coefficients in the potential. As noted previously, any change in the constants defining 
the ranges of the quadratic and cubic terms can be compensated by a rescaling of A and v. 
The dependence on the quartic coefficient range is illustrated in Fig. [TJ where we plot the 




10 



c/1 








'o 




Ph 






100 














w 
'■^ 


30 


ITS 






10 






o 






3 






1 


1 







123456789 10 
N 



FIG. 3: The median number of saddle points for quartic potentials, with the bars indicating the 
range from the 25th to the 75th percentile. 
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FIG. 4: The same as in Fig. [U but for cubic potentials. 
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FIG. 5: The same as in Fig. [2l but for cubic potentials. 
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FIG. 6: The same as in Fig. [3l but for cubic potentials. 
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FIG. 7: The dependence of the median value of s on the range parameter 04. The data shown are 
for N = 2. 

median value of s as a function of 04 with N = 2. We see that increasing 04 produces a 
roughly exponential decrease in the median value of s; decreasing 04 to zero simply reduces 
to the purely cubic potential. Hence the net effect of increasing the range for 04 would be to 
lead to a lower tunneling exponent, thus strengthening the effects that we find. The results 
for the other quantities and other values of give similar results. 

Our data indicate that the median value of s falls rapidly as the number of fields is 
increased. From this we can conclude that at large the typical local minimum is less 
likely to have a low nucleation rate and a long lifetime. This by itself does not tell us 
about the number of outliers, vacua with high values of s. To study this we examined the 
distributions of the quantities that we have plotted within a given ensemble. For example, 
in Fig. [8] we show the distributions of values of s for several choices of N with a quartic 
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potential. These roughly coincide when they are plotted as functions of s/smcdian- The data 
suggest that the frequency of large s has an approximately exponential falloff that can be 
described by 

n(s) ^ no exp (-7s/ Smedian) • (3.3) 

The values of 7 for various values of are shown in Table [Til Similar results are found for 
the data on the height of and distance to the lowest saddle point. 
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TABLE II: Best fit values for the parameter 7 in the distributions of surface tension s for cubic 
and quartic non-supersymmetric and quartic supersymmetric potentials with various A^. 

Inserting the fit of Eq. (13.21) for Smedian into this expression gives 

n{s) ^ no exp (-—^N'^s) , (3.4) 

\ tension / 

where for convenience we have defined a = atonsion- Extrapolating to larger N and using 
the estimate in Eq. (12. 9p suggests that the fraction of potentials with a tunneling exponent 
greater than some value 13 is roughly 

/(S)~exp(-/3Ar"5), (3.5) 
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FIG. 8: The distribution of values of s in the ensembles for a quartic potential with various values 
of N. The vertical axis represents the natural logarithm of the number of values in each bin. 
For the sake of clarity, the data for different values of N have been offset by constants, so only 
the slopes are meaningful. Purple diamonds correspond to = 2; green down-pointing triangles, 
= 4; red up-pointing triangles, = 6; blue squares, N = 8; and black circles, N = 10. 



where 



(3 = 10" 



7A 



(3.6) 



r' 

'-'tension 

Our numerical results suggest that 7/Ctension is close to unity. Because A is extracted from 
the shape of the landscape, there is no reason to expect it to be a small coupling constant. 
It could well be of order unity, in which case /3 ~ 10^^. Equation f l3.5p then represents 
a tremendous suppression with increasing of vacua with nucleation rates low enough to 
maintain metastability. 

Note that we have ignored gravitational effects in our analysis of tunneling rates. With 
low barriers and rapid bubble nucleation this is generally a good approximation. It should 
be noted, though, that in a de Sitter vacuum with a relatively fiat barrier the Hawking-Moss 
bounce provides an alternative mode of vacuum decay. The corresponding decay exponent 
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where Vfv and Vbarr are the values of the potential in the false vacuum and at the top of 
the barrier, while Gat is Newton's constant. Our numerical studies indicate that the barrier 
height relative to the false vacuum (i.e., Vbarr — V^) has a power law falloff with quite 
similar to that found in the effective wall tension s, thus implying a similar enhancement of 
the Hawking-Moss transition rate. 

IV. POSSIBLE IMPLICATIONS FOR THE STRING LANDSCAPE 

The results we have described so far are a general feature of quantum field theories. Our 
motivation for studying random multi-field potentials, however, is the landscape of string 
theory. In this section, therefore, we assess the applicability and implications of our findings 
for the string landscape. 

Not only does the enormous number of local minima in the high dimensional string 
theory landscape (which can serve as the endpoints of tunneling events) offer the possibility 
of significant enhancement of tunneling probabilities,^ but we have also seen that such large 
numbers of fields result in the exponential suppression of tunneling barriers. The question 
we now briefly consider is the extent to which the random potentials we have analyzed in 
the previous sections provide accurate insight into properties of the string landscape. 

For definiteness, focus on a standard flux compactification of type IIB string theory on 
a Calabi-Yau manifold M. With G = F - tH, F and H the RR and NS-NS three-form 
fluxes, and r the axio-dilaton, the Gukov-Vafa-Witten superpotential W is given by 24 1 



W = J GAn, (4.1) 
where fl is the holomorphic 3-form on M. The associated flux potential V is given by 

VM = e^{DpWDf'W -3\W\'^). (4.2) 



^ The density of vacua near the conifold point, pconifoid, in a one-dimensional moduH space is described by 
l/[r^(C + logr)^], where r is the distance from the conifold point [isf; applying this result near a generic 
point along the conifold locus in an n-dimensional moduli space shows the rapid growth in the number of 
vacua, / d'V pconifoid , with n. 
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Assume that G, together with additional contributions (e.g. D3-brane instantons or 
wrapped D7-branes), stabihzes all moduli 23|. Then consider {Va/} as M varies among 
all Calabi-Yau's, G varies over flux values, the Kahler stabilizing contributions vary over 
all possibilities, and supersymmetry-breaking effects similarly scan a broad range. With all 
such variations in play, the local minima of {Va/} will sweep through a large class of vacua. 
In the vicinity of any such vacuum state, we can expand the effective potential, yielding a 
fleld theory model of the form (12. 2p . Now, considering all such expansions that arise from all 
local minima of the collection {Vm}, we expect the expansion coefficients to randomly vary. 
Indeed, such random variation inspired the ansatz we chose in our numerical studies. But 
in applying the results of the previous sections to the string landscape, there are a number 
of details and complications that deserve further attention. 

First, we have assumed that the action for tunneling trajectories is well approximated by 
our proxy: twice that of a straight line in fleld space connecting the minima being studied 
and the optimal saddle point on the surrounding barrier. Yet, tunneling trajectories in 
multidimensional fleld spaces are notoriously subtle and can exhibit unexpected features; 
explicit examples in the landscape are the conifunneling trajectories between monodromy- 



related flux vacua found in 



25| (see also 26|, l27|]). An additional, and potentially pivotal. 



complication is that the different flux vacua that are not monodromy related are generally 
minima of distinct potentials. Physically, such transitions invoke features not captured by 
our local fleld theoretic model, including for example the nucleation of branes to absorb 



changes in flux 



28 



29| . These effects might signiflcantly affect the tunneling action, and 



possibly mitigate the fleld theory instabilities we have identifled. 

Second, we have assumed that as the stabilizing contributions to a given model are varied 
the effective potentials around local minima will have expansions that are well modeled by 
random polynomials. Is this correct? To be concrete, consider the part of the potential 
arising from the flux, G. On a given Calabi-Yau, M, there are dim H^{M) fluxes that enter 
G. As those fluxes vary, we have far fewer free parameters than the coefficients in Eq. (12. 2p . 
However, the associated minima of V will occur at different locations, p, in the moduli space 
of M. As the holomorphic form Q depends on p, the coefficients in a local expansion of the 
superpotential, Eq. (14. ip . and the corresponding potential, Eq. (14. 2p . will also vary with p. 
So, the local potential will be randomized both by the varying fluxes and by the varying 
values of the period vector over the moduli space. It seems reasonable to us that this will 
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result in local expansions well modeled by the random potentials invoked in Sec. Illlt but we 
do not have a firm argument. We will return to this issue in a forthcoming work, where we 
will explicitly check this for specific Calabi-Yau's with low- dimensional moduli spaces. 

Third, we have taken a canonical form for the kinetic terms in our field theoretic models. 
It is well known, however, that string vacua are densest in the vicinity of the conifold locus, 
where the classical moduli space metric suffers from a curvature singularity. In particular, 
near a generic point on the conifold locus we can choose local coordinates [Z^, Z^, . . . , Z^) 
on the moduli space such that Z^ = labels the conifold. Near Z^ = 0, the moduli space 
metric G behaves as: 

-IndZ^n. (4.3) 
The local form of the action then takes the form 

J V^lg'^G^j d,Z' d,Z^ - V{Z)], (4.4) 

where g is the space-time metric. In this expression the coordinates Z are the moduli space 
representation of the scalar fields (p and V is their fiux potential. At any non-singular point 
we can, of course, use a local change of field variables to absorb G into the Z, yielding 
a canonical kinetic term. But as we approach a conifold point, this change of variables 
corresponds to reducing the barrier heights in V (assuming V is continuous and is being 
expanded about a local minimum) and thus increases tunneling rates. In the regions of 
moduli space that are most densely populated with string vacua, we therefore expect the 
non-canonical kinetic terms to augment the destabilization we have found. 

Fourth, since the fiux potential V is derived from the superpotential W, and it is W 
that directly incorporates fiux values, a more accurate representation of the landscape arises 
from randomly varying coefficients in a local expansion of W, and then using the result to 
calculate the random potentials. In principle, the relationships between the coefficients in 
V, which reflect its origin in W, could alter our findings. 

We have undertaken such an analysis. Specifically, we considered a theory with chiral 
superfields $j. The superpotential W was taken to be a polynomial 

^ = lll cS^^^ + 1 E ^S*^*.*^ + \ E ^S^*^*.*^*^ + - (4-5) 

i ijk ijkl 
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The were randomly chosen real numbers in the range [0, 1], while C^jl and Cljl^ were 
complex numbers whose real and imaginary parts were taken randomly from the interval 
[— 1, 1]. The scalar field potential derived from W was then truncated at quartic order and 
analyzed in a fashion similar to our non-supersymmetric potentials. With up to = 5 
superfields we again find a power law falloff with A^. Because the N chiral superfields 
correspond to 2N real scalar fields, instead of writing the fits to the data as in Eq. f l3.2p we 
write, e.g., 

Smcdian = Ctcnsion(2iV)-". (4.6) 

As can be seen from the data in our tables, the various fit parameters are rather similar to 
those for the nonsupersymmetric case, with the most notable difference being that Smedian 
varies as A^~'^'^^, compared to A^~^-^^ in the non-supersymmetric case^. 

This example also illustrates the challenges of investigating random potentials for large 
numbers of fields. In all of our studies, supersymmetric or not, computational considerations 
have forced us to only probe a limited range of values of A^, the number of fields, and work 
with truncated potentials. We are assuming that the pattern we have found, as evidenced 
in Figs. [iHHl will continue to hold as these constraints are relaxed. 

V. A MULTIVERSE EXPLANATION OF THE COSMOLOGICAL CONSTANT? 

Our results have immediate implications for attempts to find a multiverse solution to 
the cosmological constant problem in such field theory models. The argument for such a 
solution is that even if the natural scale for the cosmological constant in a typical vacuum 
is Planckian, one might expect to find vacua with A ~ 10~^^° in Planck units if the number 
of vacua is much greater than 10^^*^. Because we want a vacuum that has not only a small 
A but also a long lifetime, what we actually need is that the number of truly metastable 
vacua be much greater than 10^^*^. 

Let us suppose that the number of vacua is A/'vac ~ J^^ , with perhaps J-" = 10. For 
metastability, we require that the tunneling exponent B be no smaller than a value -Bmm of 



^ Because we have not included any supersymmetry-breaking terms in the potential, the vacua here are 
actually stable, and our data correspond to domain walls rather than bubble walls. However, we do not 
expect the picture to be changed materially when the vacua are lifted. 
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order unity. The number of metastable vacua is then 

Ar.ac/(i?mi„)~-F^e-^^-^'^ (5.1) 

with f{B) given by Eq. (13. 5p . The requirement that this be greater than 10^^° can be written 

as 

with 

b = 13B^^ = {^^\ A5,,i, ~ 10-=^A . (5.3) 

V '-^tension / 

This clearly needs large A^, at least more than 120 if J-" = 10 as in the usual analysis. 
The new feature here is that, because of the enhancement of the tunneling rate at large N , 
taking to be too large makes matters worse rather than better. In fact, there may be no 
value of for which this condition is satisfied. This is illustrated in Fig. |9l where we have 
taken a = 2.66, the power we obtained from the analysis of quartic non-supersymmetric 
potentials. For J-" = 10 the allowed values of and b correspond to the region to the left of 
the solid curve. We see that there are no acceptable values of A^ for b > 1.4 x 10~^, and that 
the range of allowed A^ is relatively restricted until b falls well below this value. If instead 
we take J-" = 100, the allowed region extends to the dashed line, with the maximum allowed 
b increased by roughly a factor of four. 

The effect of varying the power a can be seen in Fig. [101 where we have set a = 3.16, the 
value from our supersymmetric data. Although the curve is similar to that in the previous 
case, the values of b have fallen by more than an order of magnitude. 



VI. CONCLUSIONS 



Motivated by the string landscape, we have undertaken a study of the stability of vacua in 
multi-field quantum theories. Our study has focused on random potentials with the number 
of fields A^ < 10 and on polynomial expansions that include no more than quartic field 
contributions. Even with these constraints, coming from computational feasibility, the data 
we have accumulated provide evidence that transition rates are so rapidly enhanced as a 
function of A^ that all but an exponentially small fraction of generic local minima in such 
quantum field theories are unstable to rapid decay. 
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FIG. 9: Parameter ranges allowing multiverse explanations of the cosmological constant. With 
a = 2.66 and J- = 10, a sufficient number of metastable vacua is only possible for parameters in 
the region to the left of the solid line. This region is extended to the dashed line if = 100. 




2.x 10"^ 4.x 10"^ 6.x 10"^ 8.x 10"^ 0.00001 

b 



FIG. 10: Like Fig. El but for a = 3.16 and T = 10. 
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One consequence is that the range of parameters that give a sufficient number of 
metastable vacua to provide a natural solution to the cosmological constant problem is 
severely restricted. For example, if the moduli space dimension N = 500 and there are 
J^^ = 10^°° vacua, the "coupling" A characterizing the landscape must be less than 1/20 or 
so; with = 2000, A must be an order of magnitude smaller. Alternatively, if A = 1/2 and 
N — 500, then must be greater than 6 x 10^. 

These considerations are potentially relevant to any model invoking anthropic explana- 
tions for the cosmological constant in which the required diversity of vacua is due to the 
model containing a high dimensional moduli space. However, since our results have been 
obtained in the context of a field theory with multiple scalar fields, the impact of simi- 
lar considerations on the ability of the string landscape to offer a natural solution to the 
cosmological constant puzzle will require further study. 
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